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Abstract 

SESOIL  is  a  seasonal  soil  compartment  model.  A  one-dimensional  vertical 
transport  model  is  designed  to  predict  seasonal  solute  distribution  in  the  soil 
profile  and  watershed.  CRREL  received  a  request  from  the  U.S.  Army  Engineer 
District  Alaska,  Environmental  Technical  Engineering  Office,  to  provide  techni¬ 
cal  assistance  in  evaluating  the  SESOIL  model  for  helping  to  assess  benzene 
leaching  in  the  Alaskan  environment.  The  major  tasks  outlined  in  the  request 
were  as  follows:  work  an  example  problem  for  a  diesel-contaminated  site,  do 
analytical  checks  and  do  manual  SESOIL  analytical  calculations  for  one  cycle. 
The  SESOIL  model  requires  57  input  variables  supplied  by  the  user.  An  addi¬ 
tional  8  parameters  are  required  for  the  execution  file.  This  study  did  sensitivity 
analyses  on  soil  bulk  density,  intrinsic  permeability,  disconnectedness  index, 
porosity,  organic  carbon,  adsorption  coefficient  on  organic  carbon,  and  bio¬ 
degradation  rates  of  solid  and  liquid  phases.  The  model  is  very  sensitive  to  all 
the  parameters  studied.  Despite  its  several  difficulties,  the  model  is  popular 
among  regulators  and  users  because  of  its  simplicity  compared  to  research 
models.  It  can  be  used  as  a  screening-level  tool  in  assessing  chemical  move¬ 
ment  in  the  soil  column  with  considerable  site-specific  calibrations. 


For  conversion  of  SI  units  to  non-SI  units  of  measurement  consult  Standard 
Practice  for  Use  of  the  Internafional  System  of  Units  (SI),  ASTM  Standard  E380- 
93,  published  by  the  American  Society  forTesting  and  Materials,  1916  Race  St., 
Philadelphia,  Pa.  19103. 
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Evaluating  the  SESOIL  Model  for 
Benzene  Leaching  Assessment  in  Alaska 

GURDARSHAN  S.  BRAR 


INTRODUCTION 

The  United  States  Army  Engineer  District,  Alas¬ 
ka,  Environmental  Technical  Engineering  Office, 
asked  CRREL  for  technical  assistance  in  evaluating 
the  SESOIL  model  for  aiding  in  assessing  benzene 
pollutant  leaching  in  Alaska.  The  major  tasks  out¬ 
lined  in  the  request  were  1)  to  work  an  example 
problem  for  a  diesel-contaminated  site,  2)  do  ana¬ 
lytical  checks  and  3)  make  manual  SESOIL  an¬ 
alytical  calculations  for  one  cycle. 

Drewett  et  al.  (1993)  used  the  SESOIL  model  at 
Fort  Greely  to  assess  the  presence  and  concentra¬ 
tions  of  petroleum  hydrocarbons  in  the  soil  and 
groundwater  at  five  sites.  In  this  study,  a  SESOIL 
sensitivity  analysis  was  made  for  benzene  con¬ 
tamination  in  Building  110  at  Fort  Greely  by  using 
input  data  files  of  Drewett  et  al.  (1993). 

SESOIL  history 

SESOIL  stands  for  seasonal  soil  compartment 
model.  It  is  a  one-dimensional  vertical  transport 
model  for  the  unsaturated  soil  zones  and  is  de¬ 
signed  to  predict  solute  distribution  in  the  soil 
profile  and  watershed  seasonally.  In  1981,  Arthur 
D.  Little,  Inc.  (ADL)  developed  the  program  for  the 
Environmental  Protection  Agency's  (EPA)  Office 
of  Water  and  Office  of  Toxic  Substances  (OTS).  The 
original  code  included  three  soil  columns.  In  1984, 
ADL  updated  the  model  and  included  a  fourth  soil 
layer  as  well  as  soil  erosion  algorithms  (Bonazoun- 
tas  and  Wagner  1984).  In  1985,  the  EPA's  Environ¬ 
mental  Research  Laboratory  in  Athens,  Georgia, 
thoroughly  evaluated  SESOIL  (Watson  and  Brown 
1985),  found  numerous  deficiencies  in  the  model 
and  recommended  improvements.  Furthermore, 
the  model  was  extensively  tested  at  the  Oak  Ridge 
National  Laboratory  (ORNL)  and  modified  to  en¬ 
hance  its  capabilities  (Hetrick  et  al.  1986,  1989; 


Hetrick  and  Travis  1988).  SESOIL  was  incorpo¬ 
rated  into  a  system  called  PCGEMS  (Graphical 
Exposure  Modeling  System  for  the  PC),  a  com¬ 
plete  information  management  tool  developed 
for  EPA-OTS  and  designed  to  help  users  perform 
exposure  assessments  (General  Science  Corpora¬ 
tion  1987,  1989).  General  Science  Corporation 
(1990)  further  made  improvements  to  PCGEMS 
and  incorporated  it  into  a  new  system  called 
RISKPRO  (Hetrick  et  al.  1993). 

Objective 

The  objective  of  this  report  was  to  evaluate  the 
sensitivity  of  the  SESOIL  model  for  assessing  ben¬ 
zene  leaching  in  the  Alaskan  environment. 

LITERATURE  REVIEW 

Donigian  and  Rao  (1986a)  comprehensively  re¬ 
viewed  the  model  and  concluded  that  the  chemi¬ 
cal  transport  part  of  the  code  needed  improve¬ 
ments.  The  model  did  not  handle  the  large 
variations  in  the  hydraulic  properties  of  the  soil 
profile.  Furthermore,  they  reported  that  SESOIL 
overpredicted  the  mass  of  aldicarb  movement  to 
groundwater  (Donigian  and  Rao  1986b).  Daugh¬ 
erty  (1993)  found  that  SESOIL  does  not  simulate 
the  transfer  of  contaminants  from  the  NAPL  to 
the  aqueous  phase,  but  assumes  that  all  contami¬ 
nants  are  dissolved. 

Ladwig  et  al.  (1993)  tested  the  SESOIL  model 
for  benzene,  ethylbenzene,  toluene,  xylene  and 
1,2-dichloroethane  under  Wisconsin  conditions 
and  discovered  that  the  model  was  most  sensitive 
to  soil  type,  biodegradation  rate,  residual  concen¬ 
tration,  residual  layer  thickness  and  organic  car¬ 
bon  content  of  the  soil.  The  model  was  least  sensi¬ 
tive  to  climate,  soil  disconnectedness  index,  pH 


and  depth  of  contaminant  burial.  Furthermore, 
the  model  was  highly  sensitive  to  intrinsic  per¬ 
meability  for  benzene  in  sand.  However,  it  was 
insensitive  to  intrinsic  permeability  for  treatment 
of  benzene  in  till.  The  till  was  a  finer  grained  envi¬ 
ronment,  with  fewer  available  pore  spaces  than 
the  sand. 

Bonazountas  and  Kallidromitou  (1993)  report¬ 
ed  that  at  the  present  state  of  scientific  knowl¬ 
edge,  SESOIL  is  a  well-developed  soil  compart¬ 
ment  model.  They  emphasized  that  the  model 
must  be  calibrated  for  site-specific  variables  such 
as  biodegradation  rate,  soil  organic  carbon  con¬ 
tent,  local  climate  and  depth  to  the  water  table. 
Hetrick  et  al.  (1993)  recommended  that  predic¬ 
tions  for  the  hydrology  at  a  given  site  be  calibrat¬ 
ed  to  agree  with  known  measurements.  They 
suggested  that  the  user  must  conduct  sensitivity 
analyses  or  evaluate  results  obtained  by  assign¬ 
ing  distributions  to  the  input  parameters  (O'Neill 
et  al.  1982,  Gardner  1984,  Hetrick  et  al.  1991). 
When  properly  used,  SESOIL  is  an  effective 
screening-level  tool  for  assessing  chemical  move¬ 
ment  in  soil. 

Oregon  State  regulator  Anderson  (1992)  tested 
the  SESOIL  model  to  simulate  transport  of  organ¬ 
ic  contaminants  through  the  vadose  zone.  He  also 
used  the  MINTEQA1  model  (USEPA  1987)  to 
simulate  inorganic  speciation  and  concentrations 
under  specified  field  conditions.  Sensitivity  anal¬ 
yses  were  performed  on  organic  contaminants — 
benzene,  carbon  tetrachloride,  ethylbenzene, 
methyl  chloride,  naphthalene,  perchloroethylene, 
toluene,  1,1,1-trichloroethane,  1,1,2-trichloro- 
ethylene  and  o-xylene.  The  results  of  the  tests  in¬ 
dicated  that  the  most  important  parameters  con¬ 
trolling  organic  contaminant  transport  in  the 
SESOIL/ AT123D  models  were  the  thickness  of 
the  contaminated  zone,  depth  from  the  contami¬ 
nated  zone  to  groundwater,  hydraulic  conductiv¬ 
ity  and  gradient,  and  the  fraction  of  the  organic 
carbon  in  the  soil. 

Numerous  studies  have  been  conducted  to 
compare  the  sensitivity  of  the  SESOIL  model  with 
other  models,  using  field  input  data  (Bonazoun¬ 
tas  et  al.  1982;  Wagner  et  al.  1983;  Hetrick  1984; 
Kincaid  et  al.  1984;  Watson  and  Brown  1985;  Mel- 
ancon  et  al.  1986;  Hetrick  et  al.  1986, 1989;  Hetrick 
and  Travis  1988).  The  model  was  used  in  risk 
assessments  to  evaluate  direct  coal  liquefaction 
(Walsh  et  al.  1984),  incineration  of  hazardous 
waste  (Holton  et  al.  1985,  Travis  et  al.  1986),  trans¬ 
port  of  benzene  to  the  groundwater  (Tucker  et  al. 
1986),  soil  cleanup  levels  in  California  (Oden- 


crantz  et  al.  1991,  1992)  and  site  sensitivity  rank¬ 
ing  for  Wisconsin  soils  (Ladwig  et  al.  1993). 

Several  computer  models  are  in  use  to  assess 
the  cleanup  levels  of  hydrocarbon  contaminated 
soils.  Recently,  Odencrantz  et  al.  (1992)  reviewed 
the  following  vadose  zone  transport  models: 

1.  PRZM — Pesticide  Root  Zone  Model 

2.  CMLS — Chemical  Movement  in  Layered 
Soils 

3.  GLEAMS — Groundwater  Loading  Effects  of 
Agricultural  Management  Systems 

4.  LE  ACHMP — Leaching  Estimation  And  Chem¬ 
istry  Model-Pesticide 

5.  MOUSE — Method  Of  Underground  Solute 
Evaluation 

6.  PESTAN — Pesticide  Analytical  model 

7.  Jury's  analytical  transport  model 

8.  MOFAT — Multiphase  Organic  Flow  And 
Transport 

9.  SESOIL. 

They  recommended  the  SESOIL  and  MOFAT 
models  on  the  basis  of  the  following  vadose  zone 
transport  criteria: 

1.  Needed  is  a  one-dimensional  vertical  model 
that  can  be  discretized  into  at  least  three  lay¬ 
ers,  with  varying  soil  properties  and  chemi¬ 
cal  loading  concentrations  by  layer. 

2.  The  model  must  represent  one  mobile  chemi¬ 
cal  component  in  three  phases  (adsorbed, 
aqueous  and  gaseous). 

3.  The  model  must  represent  biodegradation, 
volatilization  and  gaseous  diffusion. 

4.  The  model  must  represent  the  variable  effects 
of  local  climate  (precipitation,  temperature, 
etc.). 

MOFAT  is  a  two-dimensional  finite-element 
modeling  code  that  transports  up  to  five  compo¬ 
nents  between  four  phases  (air,  water,  oil  and 
soil)  and  also  allows  up  to  ten  soil  layers  of  differ¬ 
ent  properties.  MOFAT  is  more  powerful  but 
harder  to  use  than  SESOIL;  therefore,  fewer  users 
have  adopted  it. 

The  Council  for  Health  and  Environmental 
Safety  of  Soils'  (CHESS)  Analysis  and  Environ¬ 
mental  Fate  Committee  reviewed  several  existing 
environmental  fate  models  and  recommended 
the  SESOIL  and  POSSM  models  for  further  con¬ 
sideration  as  tools  to  aid  the  risk  assessment  pro¬ 
cess  (Calabrese  and  Kostecki  1992).  POSSM, 
which  stands  for  PCB  On-Site  Spill  Model,  is  a 
contaminant  transport  model  developed  to  pre¬ 
dict  environmental  concentrations  associated 
with  a  chemical  spill.  It  is  a  modified  version  of 
the  PRZM  model  and  is  considered  as  a  one- 
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dimensional,  compartmental,  dynamic  transport 
and  fate  model.  The  model  simulates  daily 
changes  in  chemical  concentrations  on  a  spill 
site's  soil  and  vegetation,  as  well  as  losses  of 
chemicals  attributable  to  volatilization,  surface 
runoff-soil  erosion  and  leaching  to  groundwater. 

Vertical  And  Lateral  Organic  Redistribution 
(VALOR),  a  two-dimensional  numerical  model 
for  immiscible  multiphase  fluid  flow  in  subsur¬ 
face  systems,  was  developed  by  Abriola  et  al. 
(1992).  VALOR  solves  only  flow  (pressure)  equa¬ 
tions  for  different  phases  of  air,  water  and  organic 
liquids;  however,  MOFAT  solves  both  flow  and 
mass  transport  (concentration)  equations. 

SESOIL  EVALUATION 
SESOIL  processes 

SESOIL  simulates  the  major  processes  that  can 
affect  pollutant  transport  in  the  unsaturated 
zone,  including  moisture  movement,  volatiliza¬ 
tion  and  vapor  phase  diffusion,  adsorption,  and 
chemical  as  well  as  biological  degradation  (Lad- 
wig  et  al.  1993).  The  hydrological  cycle  compo¬ 
nent  is  based  on  the  theory  of  Eagelson  (1978)  that 
uses  a  statistical  approach  to  water  balance  calcu¬ 
lations  for  estimating  one-dimensional  soil-water 
movement.  Parameters  considered  in  the  hydro- 
logic  cycle  are  rainfall,  surface  runoff,  infiltration, 
soil-water  content,  evapotranspiration  and 
groundwater  runoff.  The  sediment  cycle  relates 
to  sediment  washload  as  a  result  of  rainstorms 
(i.e.,  soil  erosion  from  surface  runoff).  To  incorpo¬ 
rate  the  effect  of  washload,  the  sediment  module 
uses  a  combined  statistical  and  deterministic  ap¬ 
proach.  The  pollutant  fate  cycle  includes  convec¬ 
tive  transport,  volatilization,  adsorption-desorp¬ 
tion,  chemical  degradation-decay,  biological 
transformation,  hydrolysis,  photolysis,  oxidation 
and  complexation  with  metals.  This  module  uses 
a  deterministic  approach  with  analytical  equa¬ 
tions  describing  chemical  and  biological  process¬ 
es  that  affect  contaminant  transport.  Details  of 
these  cycles,  along  with  the  transport  equations, 
are  given  in  the  New  SESOIL  User's  Guide  (Hatrick 
et  al.  1993). 

SESOIL  assumptions 

The  SESOIL  model  is  based  on  three  assump¬ 
tions: 

1.  Pollutant  concentrations  in  all  phases  and  in 
all  compartments  of  the  soil  system  are  at 
equilibrium  at  all  times. 

2.  The  law  of  mass  conservation  determines  the 


equilibrium  concentration  of  chemical  spe¬ 
cies  over  a  series  of  monthly  or  yearly  steps. 

3.  Pollutant  transport  takes  place  in  the  unsatu¬ 
rated  soil  zone. 

SESOIL  input  parameters 

SESOIL  has  options  for  annual  or  monthly  sim¬ 
ulations.  The  annual  simulations  need  annual  cli¬ 
matic  data  and  monthly  simulations  require 
monthly  input  data.  The  annual  simulation  op¬ 
tion  is  not  available  in  the  new  RISKPRO  system; 
therefore,  it  is  not  included  in  this  report  with  the 
exception  of  the  hydrological  cycle,  which  imple¬ 
ments  the  annual  algorithm.  The  monthly  input 
data  are  grouped  into  four  types:  climate  data, 
soil  data,  chemical  data  and  initial  chemical  dis¬ 
tribution  in  the  soil.  The  user  must  supply  57  in¬ 
put  variables  to  run  the  model  (Ladwig  et  al. 
1993):  12  monthly  input  values  are  required  for  33 
of  the  parameters.  In  the  execution  file,  8  addi¬ 
tional  parameters  are  required.  The  washload  in¬ 
put  file  was  not  used  at  the  Fort  Greely  site;  there¬ 
fore,  a  discussion  of  the  sediment  module  is  not 
included.  Details  of  input  variables  for  the  four 
categories  are  as  follows. 

Climate  data 

1.  Average  monthly  temperature. 

2.  Average  monthly  fraction  of  cloud  cover. 

3.  Average  monthly  relative  humidity. 

4.  Average  monthly  shortwave  albedo. 

5.  Monthly  precipitation  depth. 

6.  Mean  storm  duration. 

7.  Number  of  storms. 

8.  Latitude. 

9.  Mean  length  of  rain  period. 

10.  Average  daily  evapotranspiration  rate. 

Soil  data 

1.  Description  of  soil  type. 

2.  Soil  bulk  density. 

3.  Intrinsic  permeability. 

4.  Disconnectedness  index. 

5.  Effective  soil  porosity. 

6.  Organic  carbon  content. 

7.  Freundlich  exponent. 

8.  Cation  exchange  capacity. 

Chemical  data 

1.  Description  of  the  compound. 

2.  Solubility  in  water. 

3.  Diffusion  coefficient  in  air. 

4.  Henry's  Law  constant. 

5.  Adsorption  coefficient  on  organic  carbon. 

6.  Overall  adsorption  coefficient. 

7.  Molecular  weight  of  the  compound. 
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8.  Valence  of  the  compound. 

9.  Neutral  hydrolysis  constant. 

10.  Base  catalyzed  hydrolysis  constant. 

11.  Acid  catalyzed  hydrolysis  constant. 

12.  Biodegradation  rate  in  the  liquid  phase. 

13.  Biodegradation  rate  in  the  solid  phase. 

14.  Stability  constant  of  compound-ligand  com¬ 
plex. 

15.  Number  of  moles  of  ligand  per  mole  of  com¬ 
pound  complexed. 

16.  Molecular  weight  of  ligand. 

Application  data 

1.  Number  of  soil  layers. 

2.  Number  of  years  of  data  included  in  the  data 
set. 

3.  Surface  area  of  the  compartment. 

4.  Depth  of  each  layer. 

5.  pH  of  each  layer. 

6.  Intrinsic  permeability  of  each  layer. 

7.  Layer  ratios  for  biodegradation  in  the  liquid 
phase. 

8.  Layer  ratios  for  biodegradation  in  the  solid 
phase. 

9.  Layer  ratios  for  organic  carbon  content. 

10.  Layer  ratios  for  cation  exchange  capacity. 

11.  Layer  ratios  for  the  Freundlich  exponent. 

12.  Layer  ratios  for  adsorption  coefficient. 

13.  Monthly  pollutant  load  (mass/unit  area) 
entering  each  zone. 

14.  Monthly  mass  of  pollutant  transformed  in 
each  sublayer  by  some  other  process. 

15.  Monthly  mass  of  pollutant  removed  from 
each  sublayer  by  some  other  process 

16.  Monthly  ligand  mass  input  to  each  sublayer. 

17.  Index  of  volatilization-diffusion  occurrence 
from  each  layer. 

18.  Index  of  subsurface  pollutant  runoff. 
RISKPRO  contains  a  file  management  system  that 

simplifies  data  input  for  SESOIL.  Default  climatic 
data  for  a  particular  state  can  be  purchased  and  ac¬ 
cessed  through  the  RISKPRO  package. 

SESOIL's  merits 

Being  a  part  of  PCGEMS  and  now  RISKPRO,  SE¬ 
SOIL  is  popular  among  regulators  for  risk  as¬ 
sessment  studies.  The  compartmental  module  lets 
many  users  run  the  model  for  specific  data  sets  or 
site  conditions.  The  model  had  continuous  sup¬ 
port  from  EPA-OTS.  Periodically,  SESOIL  is  im¬ 
proved  and  modified  by  the  Oak  Ridge  National 
Laboratory. 

The  scientific  community  engaged  in  chemical 
fate  modeling  has  accepted  and  recognized  the 


SESOIL  model.  Extensive  validations  that  work 
under  different  site  specific  scenarios  must  be 
done  to  calibrate  the  model.  The  model  is  struc¬ 
tured  to  simulate  chemical  transport  for  more  than 
a  month.  SESOIL  accommodates  physical,  chemi¬ 
cal  and  biological  changes  in  the  contaminant 
leaked  into  the  soil  system. 

The  sensitivity  analyses  were  conducted  on  ad¬ 
sorption  and  volatilization  for  different  soil  types 
and  in  different  climates.  SESOIL's  hydrologic  cy¬ 
cle  has  been  found  to  be  a  good  long-term  predic¬ 
tor  for  groundwater  and  surface  runoff,  evapo- 
transpiration  and  infiltration.  Uncertainty  analysis 
is  introduced  into  the  hydrological  cycle  with 
probability  density  functions,  which  produces 
probability  distributions  of  water  balance  and 
yields  long-term  seasonal  averages  of  the  water 
balance. 

SESOIL's  demerits 

Application  of  the  SESOIL  model  in  the  Alaskan 
environment  requires  extensive  calibration  and 
modifications  because  it  does  not  function  at  tem¬ 
peratures  below  freezing  (Calabrese  and  Kostecki 
1992).  The  present  code  for  SESOIL  uses  a  single 
homogeneous  soil  column  for  the  hydrological  cy¬ 
cle.  In  the  Alaskan  environment,  because  soils  usu¬ 
ally  have  discontinuous  and  fractured  permafrost 
underneath,  the  model  will  need  modifications. 
SESOIL  will  not  work  at  sites  having  large  vertical 
variations  in  soil  properties. 

The  use  of  SESOIL  is  limited  because  it  requires 
site-specific  data  for  calibration.  When  input  data 
are  not  available,  the  model  user  must  use  com¬ 
plex  calculations  to  generate  an  input  file.  The  user 
needs  the  expertise  to  select  the  appropriate  equa¬ 
tions  required  for  developing  input  data.  In  a  situ¬ 
ation  where  site-specific  data  are  not  available,  us¬ 
ing  data  published  in  the  literature  or  default 
values  might  simulate  results  that  are  inaccurate 
by  orders  of  magnitude. 

The  model  does  not  address  the  free  product 
movement  in  the  vadose  zone,  such  as  depth  of 
penetration  of  bulk  hydrocarbons,  spread  and  mi¬ 
gration  rate  of  free  product,  effects  of  large  con¬ 
centrations  of  other  organics  on  adsorption  and 
mobility,  and  emission  rate  for  a  pure  bulk  hydro¬ 
carbon  on  a  soil  surface.  Furthermore,  the  model 
accommodates  the  migration  of  a  single  solute  in 
an  aqueous  phase  rather  than  a  nonaqueous 
phase.  The  model's  inability  to  distinguish  be¬ 
tween  the  NAPL  phase  or  the  water  phase  (domi¬ 
nant  transport  carrier)  can  create  significant  errors 
in  simulations. 
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Simulations  with  SESOIL  are  very  poor  when 
run  with  laboratory-generated  input  data.  For  in¬ 
stance,  with  laboratory  input  data,  the  model 
simulated  that  65.5%  of  the  chemical  released  mi¬ 
grated  to  the  groundwater;  conversely,  field  data 
showed  only  4  to  8%  of  the  chemical  migrating  to 
the  groundwater  (Calabrese  and  Kostecki  1992). 
Since  these  studies  were  conducted,  the  model 
has  been  improved  and  tested  in  laboratory  col¬ 
umns  with  six  organic  compounds,  and  results 
were  compared  with  three  field  studies.  The 
modified  model  delivered  better  simulation  re¬ 
sults  for  some  compounds,  but  there  was  no  im¬ 
provement  for  compounds  with  the  lowest  and 
highest  adsorption  coefficients. 

The  current  version  of  the  SESOIL  model  can 
simulate  vadose  zone  contaminant  concentration 
in  four  soil  columns,  with  each  column  divided 
into  ten  sub-layers.  SESOIL's  hydrological  cycle 
considers  the  soil  column  as  one  homogeneous 
compartment.  Therefore,  only  one  set  of  soil- 
water  content,  porosity  and  core-disconnected¬ 
ness  parameters  is  used  to  describe  the  entire 
vadose  zone.  The  model  will  not  work  if  the  soil 
has  low  permeability.  SESOIL  does  not  address 
the  soil-water  spatial  variability  and  water  flow 
in  each  compartment.  Therefore,  chemical  trans¬ 
port  and  distribution  in  the  soil  column  could  be 
affected  because  the  retardation  coefficients  and 
volatilization  fluxes  of  certain  chemicals  are 
water-sensitive. 

The  model  considers  the  internal  soil  moisture 
at  the  beginning  of  each  storm  and  in  the  inter¬ 
storm  periods  to  be  uniform,  at  its  long-term 
space-time  average.  This  assumption  may  be  a 
considerable  departure  from  reality,  as  the  soil 
moisture  profile  in  a  later  modeling  period  (i.e., 
1  month)  cannot  be  influenced  by  the  soil  mois¬ 


ture  profile  resulting  from  a  preceding  modeling 
period.  Furthermore,  the  effects  of  snow  and  ice 
melt  on  hydrological  and  sediment  washload  cy¬ 
cles  under  Alaskan  conditions  need  to  studied. 

Validity  of  the  sediment  cycle  has  not  been  as 
thoroughly  tested  as  other  cycles  have.  This  cycle 
ignores  the  chemical  transport  while  simulating 
sediment  transport.  The  major  limitation  of  the  lay¬ 
ered  approach  is  that  when  a  chemical  enters  a  soil 
layer,  SESOIL  considers  it  to  be  uniformly  distribut¬ 
ed  throughout  the  layer.  For  instance,  the  model 
simulates  low  chemical  concentration  for  large  soil 
layers  and  high  concentration  for  small  soil  layers. 

Hetrick  et  al.  (1989)  reported  that  modified 
SESOIL  predictions  are  in  good  agreement  with  ob¬ 
served  laboratory  and  field  data.  However,  the 
model  underestimates  the  concentration  near  the 
soil  surface  and  they  speculated  that  this  may  be  a 
result  of  SESOIL  ignoring  the  upward  movement  of 
the  chemical  with  the  upward  movement  of  water 
ascribable  to  soil  evaporation  losses.  The  model  ig¬ 
nores  the  diffusive  mobility  of  chemicals,  which 
may  be  important,  depending  on  the  Henry's  Law 
constant  of  the  compound,  and  does  not  consider 
the  volatilization  enhancement  when  water  evapo¬ 
rates.  SESOIL  should  not  be  applied  to  a  specific 
site  with  only  limited  calibrations  (Hetrick  et  al. 
1989). 

SENSITIVITY  ANALYSIS 

The  preliminary  sensitivity  analysis  of  SESOIL 
for  benzene  was  done  by  varying  the  few  soil  and 
chemical  input  parameter  values  (Table  1)  to  below 
or  above  the  original  values  used  by  Drewett  et  al. 
(1993).  The  input  data  file  (SSOUT034.OUT)  devel¬ 
oped  by  Drewett  et  al.  (1993)  in  their  SESOIL  tests 
for  Fort  Greely,  Alaska,  was  used  in  the  sensitivity 


Table  1.  Summary  of  sensitivity  analysis  in  the  form  of  time  to 
predict  peak  and  maximum  leachate  concentration. 


Parameter 

Test* 

value 

Changed 

value 

Change  in  peak 
Cone .  Time 
( uv/mL )  (years) 

Soil  density  (g/cm3) 

1.92 

1.62 

+0.12 

-5 

Intrinsic  permeability  (cm2) 

lxlO-7 

lxl(H 

+1.84 

-14 

Disconnectedness  index  (unitless) 

10 

6.3 

-0.01 

+4 

Porosity  (fraction) 

0.25 

0.35 

-0.04 

+12 

Organic  carbon  content  (%) 

0.13 

0.09 

+0.056 

-10 

Adsorption  coefficient  on 

83 

69 

-0.013 

-5 

organic  carbon  (koc) 

Biodegradation  rate  in  solid 

5.48xl(H 

2.74X10-4 

+1.30 

+2 

phase  (mg/kg  per  day) 

Biodegradation  rate  in  liquid 

5.48X10-4 

2.74xlCH 

+0.28 

+1 

phase  (mg/kg  per  day) 

*  Drewett  et  al.  (1993) 
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Time  (yrs) 

a.  Using  Harding  Lawson  Associates  output  file 
SSOUT034.OUT. 


Time  (yrs) 


b.  Lower  soil  bulk  density  input  parameter  (1.62  com¬ 
pared  to  1.92  g/cm3). 


Time  (yrs) 

c.  Higher  intrinsic  permeability  input  parameter 
(1x1  Cr4  compared  to  1x1 0~7  cm2). 

Figure  1.  SESOIL  model  simulations  for 

tests.  The  sensitivity  of  the  parameters  is  de¬ 
scribed  in  terms  of  change  in  dissolved  maxi¬ 
mum  concentration  of  benzene  and  time  taken  to 
reach  the  maximum  benzene  concentration  at 
5330  cm  depth  compared  with  original  values  (Fig. 
la).  The  changes  in  benzene  concentration  as  ad¬ 
sorbed  and  as  air-in-pores  were  similar  to  dis¬ 
solved  concentration. 

Soil  bulk  density 

Soil  bulk  density  had  a  significant  effect  on  the 


Time  (yrs) 

d.  Decreased  disconnectedness  index  input  parameter 
(6.3  compared  to  10). 
benzene  concentration  vs.  time  at  5330-cm  depth. 

peak  concentration  of  benzene.  The  benzene  con¬ 
centration  increased  (0.12  pg/mL)  compared  to 
the  Drewett  et  al.  (1993)  estimates  when  the  soil 
bulk  density  was  lowered  from  1.92  to  1.62  g/ cm3 
(Fig.  lb).  Furthermore,  benzene  was  released  5 
years  earlier  into  the  groundwater  compared 
with  the  Drewett  et  al.  (1993)  estimates. 

Intrinsic  permeability 

SESOIL  is  highly  sensitive  to  changes  in  intrin¬ 
sic  permeability  from  1  x  10"7  to  1  x  10-4  cm2  (Ta- 
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Benzene  Concentration  Benzene  Concentration 


0.25).  (0.09  compared  to  0.13%). 


g.  Lower  adsorption  coefficient  of  organic  carbon  in¬ 
put  parameter  (69  compared  to  83). 


Figure  1 

ble  1,  Fig.  lc)  and  benzene  concentration  elevated 
from  0.04  to  1.84  pg/mL.  Benzene  leached  to  the 
groundwater  14  years  earlier  than  in  the  Drewett 
et  al.  (1993)  estimates.  These  results  agree  with 
the  findings  of  Ladwig  et  al.  (1993). 

Disconnectedness  index 

Decreasing  the  value  of  the  disconnectedness 
index  (10  to  6.3)  lowered  benzene  concentration 
(-0.01  pg/mL)  at  the  bottom  of  the  fourth  soil  col¬ 
umn  (Fig.  Id).  Compared  with  the  Drewett  et  al. 


h.  Lower  biodegradation  rate  in  solid  phase  input  pa¬ 
rameter  ( 2.74xl0~ 4  compared  to  5.48x10 ~4). 

(cont'd). 


(1993)  results,  benzene  travel  time  to  the  ground- 
water  is  increased  by  4  years.  These  results  are 
consistent  with  the  observations  of  Odencrantz  et 
al.  (1991)  that  SESOIL  was  highly  sensitive  to  the 
disconnectedness  index. 

Effective  porosity 

Changing  the  effective  porosity  from  0.25  to 
0.35  decreased  the  peak  benzene  concentration 
(-0.037  pg/mL)  and  increased  the  peak  leach 
time  by  12  years  (Fig.  le)  compared  to  the 
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Time  (yrs) 

i.  Lower  biodegradation  rate  in  liquid  phase  input 
parameter  (2.7ixlO~^  compared  to  SASxICt4). 

Figure  1  (coni' d).  SESOIL  model  simulations  for  ben¬ 
zene  concentration  vs.  time  at  5330-cm  depth. 

Drewett  et  al.  (1993)  estimates.  This  may  be  attrib¬ 
utable  to  the  higher  volatilization  of  the  benzene 
compound  because  of  increased  porosity,  which 
allowed  more  air-filled  voids. 

Soil  organic  carbon  content  (Foe,  %)  and  koc 
SESOIL  showed  high  sensitivity  to  organic  car¬ 
bon  content.  Organic  carbon  had  a  significant  effect 
on  the  rate  of  chemical  movement,  as  well  as  on  the 
peak  leached  concentration.  Slight  reduction  in  or¬ 
ganic  carbon  (0.13  to  0.09%)  elevated  the  benzene 
concentration  (0.56  jxg/mL,  Fig.  If),  and  the  con¬ 
taminant  was  released  to  groundwater  10  years 
earlier  than  the  original  estimates  calculated. 

The  koc,  adsorption  coefficient  on  organic  car¬ 
bon,  had  a  slight  effect  on  the  sensitivity  of 
SESOIL.  By  lowering  the  koc  value  from  83  to  69 
(Fig.  lg),  the  peak  concentration  of  benzene  ele¬ 
vated  (0.13  |ig/mL)  compared  to  the  Drewett  et  al. 
(1993)  tests.  Similarly,  koc  also  lowered  benzene 
transport  time  (-5  years)  to  the  groundwater  com¬ 
pared  with  the  original  simulations. 

Biodegradation  rate  in  solid  and  liquid  phases 

As  expected,  reducing  the  biodegradation  rate 
of  the  solid  phase  (5.48  x  10-4  to  2.74  x  10-4) 
increased  both  the  benzene  concentration  (+1.29 
|ig/mL)  and  the  travel  time  to  groundwater  (+2 
years.  Fig.  lh).  When  the  biodegradation  rates 
of  the  liquid  phase  were  reduced  to  2.74  x  10"4, 


the  benzene  concentration  was  increased  by  0.28 
pg/mL  (Fig.  li).  The  benzene  travel  time  through 
the  vadose  zone  was  also  increased  (+1  years). 

MANUAL  CHECKS  OF  OUTPUT  FILE 

For  confirmation,  the  output  results  were  man¬ 
ually  checked  for  water  balance,  benzene  mass  in¬ 
put,  benzene  mass  distribution  for  each  layer,  ben¬ 
zene  concentration  in  the  soil  water  at  the  lowest 
soil  sublayer,  and  maximum  benzene  depth.  SE¬ 
SOIL  simulated  these  parameters  precisely  as  per 
set  equation  codes. 

Water  balance 

The  water  balance  was  checked  by  using  the 
equations 


P-E- MR =S+G=Y 

(i) 

II 

i 

(2) 

where  P  =  precipitation 

E  =  evapotranspiration 
MR  =  moisture  retention 
S  =  surface  runoff 

G  =  groundwater  runoff  or  recharge 
Y  =  yield 
I  =  infiltration. 

SESOIL  provides  solutions  to  eq  1  and  2  by  iter¬ 
ating  on  soil  moisture  until  the  calculated  value 
for  precipitation  is  within  1%  of  the  user  provided 
input  value.  After  the  completion  of  this  iteration, 
the  individual  components  such  as  infiltration, 
evapotranspiration,  moisture  retention,  surface 
runoff  and  the  groundwater  runoff  (recharge)  are 
known.  Infiltration  is  checked  as  the  difference  be¬ 
tween  the  precipitation  and  the  surface  runoff, 
which  is  equal  to  the  moisture  retention  plus  the 
evapotranspiration  plus  the  groundwater  runoff. 

Benzene  mass  input 

Benzene  mass  (jag)  includes  the  amount  of  ben¬ 
zene  in  the  precipitation  and  the  amount  loaded 
into  each  of  up  to  four  major  layers.  SESOIL  com¬ 
putes  the  amount  of  benzene  in  the  precipitation 
as 

[benzene  concentration  in  rain/benzene 
maximum  solubility  in  water]  x  solubility  in 
water  x  infiltration  rate  x  area  of  application. 

The  benzene  load  in  each  layer  is  calculated  from 
the  area  of  application  x  benzene  application.  To 
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calculate  the  load  for  the  major  layer,  contained 
sublayers  are  added  to  the  first  sublayer  of  that 
particular  major  layer. 

Distribution  of  benzene  mass 

Benzene  mass  distribution  in  the  column  was 
calculated  by  adding  the  masses  of  the  contami¬ 
nant  from  all  the  phases  for  the  year  of  contami¬ 
nant  loading.  Volatilized  mass  was  added  to  the 
mass  in  moisture,  adsorbed  on  soil  and  in  soil  air 
for  each  sublayer  down  through  the  soil  column. 
The  computed  tmass  was  added  to  the  sublayer 
above  it  in  the  next  time  step.  Diffused  mass  was 
added  and  divided  by  the  time  cycle. 

Benzene  concentration  distribution 
for  lowest  sublayer 

Benzene  concentrations  (jig/mL)  were  checked 
from  the  output  file  for  each  year  until  the  benzene 
reached  the  lowest  soil  zone.  Also  recorded  was  the 
greatest  benzene  concentration  as  well  as  its  travel 
time. 

Maximum  pollutant  depth 

Drewett  et  al.  (1993)  assumed  groundwater  depth 
at  5330  cm;  therefore,  benzene  transport  was  simu¬ 
lated  at  5330  cm  as  illustrated  in  Figure  1. 

SUGGESTED  IMPROVEMENTS 

Several  workers  (Melancon  et  al.  1986;  Ander¬ 
son  1992;  Odencrantz  et  al.  1992;  Hetrick  et  al.  1989, 
1993;  Ladwig  et  al.  1993)  performed  sensitivity  anal¬ 
yses  on  SESOIL,  compared  it  with  other  vadose  zone 
chemical  transport  models,  and  modified  input  pa¬ 
rameters  according  to  the  local  conditions.  They 
emphasized  that  vadose  zone  models  need  to  be 
calibrated  with  measured  site-specific  data.  The 
suggested  improvements  for  SESOIL  calibration 
and  validation  under  Alaskan  conditions  are: 

1.  Adjust  the  input  parameters  of  SESOIL  accord¬ 
ing  to  different  climatic  regions  of  Alaska. 

2.  Propose  soil  input  parameters  as  per  the  soil 
types  of  Alaska. 

3.  Study  the  chemical  characteristics  of  different 
chemicals  spilled  at  DOD  installations. 

4.  Compare  the  SESOIL  chemical  transport  simu¬ 
lations  with  other  vadose  zone  models  avail¬ 
able  in  the  literature. 

CONCLUSIONS 

Several  computer  models  are  in  use  to  assess  the 
cleanup  levels  of  hydrocarbon-contaminated  soils. 
Vadose  zone  models  were  extensively  reviewed  by 


several  scientists  and  regulatory  agencies.  At  the 
present  state  of  scientific  knowledge,  SESOIL  is  a 
well-developed  seasonal  soil  compartment  model. 
Despite  its  several  difficulties,  the  model  is  popu¬ 
lar  among  regulators  and  users  because  of  its  sim¬ 
plicity  compared  to  research  models.  It  can  be  used 
as  a  preliminary  leachate  assessment  screening 
tool  with  considerable  site-specific  calibrations. 
Different  authors  found  different  input  parame¬ 
ters  sensitive  to  their  local  conditions  in  reference 
to  a  particular  contaminant.  The  MOFAT,  POSSM 
and  VALOR  models  can  also  be  considered  for 
leachate  assessment  in  Alaska. 
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